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ABSTRACT :In this paper, we established the condition of the occurrence of local bifurcation (such as saddle- 
node, transcritical and pitchfork)with particular emphasis on the hopf bifurcation near of the positive 
equilibrium point of ecological mathematical model consisting of prey-predator model involving prey refuge 
with two different function response are established. After the study and analysis, of the observed incidence 
transcritical bifurcation near equilibrium point E 0 , E 1 , E 2 as well as the occurrence of saddle-node bifurcation 
at equilibrium point E3.lt is worth mentioning, there are no possibility occurrence of the pitchfork bifurcation 
at each point. Finally, some numerical simulation are used to illustration the occurrence of local bifurcation o 
this model. 
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I. INTRODUCTION 

A bifurcation mean there is a change in the 
stability of equilibria of the system at the value. 
Many, if not most differential equations depend of 
parameters. Depending on the value of these 
parameters, the qualitative behavior of systems 
solution can be quite different [1]. 

Bifurcation theory studies the qualitative 
changes in the phase portrait, for example, the 
appearance and disappearance of equilibria, periodic 
orbits, or more complicated features such as strange 
attractors. The methods and results of bifurcation 
theory are fundamental to an understanding of 
nonlinear dynamical systems. The bifurcationis 
divided into two principal classes: local bifurcations 
and global bifurcations. Local bifurcations, which 
can be analyzed entirely through changes in the 
local stability properties of equilibria, periodic orbit 
or other invariant sets as parameters cross through 
critical thresholds such as saddle node, transcritical, 
pitchfork, period-doubling (flip), Hopf and Neimark 
(secondary Hopf) bifurcation. Global bifurcations 
occur when larger invariant sets, such as periodic 
orbits, collide with equilibria. This causes changes 
in the topology of the trajectories in phase space 
which cannot be confined to a small neighborhood, 
as is the case with local bifurcations. In fact, the 
changes in topology extend out of an arbitrarily 
large distance (hence “global”). Such as homoclinic 
in which a limit cycle collides with a saddle point, 
and heteroclinic bifurcation in which a limit cycle 
collides with two or more saddle points, and 
infinite-periodic bifurcation in which a stable node 
and saddle point simultaneously occur on a limit 
cycle [2]. 

the name "bifurcation" in 1885 in the first paper in 
mathematics showing such a behavior also later 
named various types of stationary points and 
classified them. Perko L. [3] established the 
conditions of the occurrence of local bifurcation 
(such as saddle-node, transcritical and pitchfork). 
However,the necessary condition for the occurrence 


of the Hopf bifurcation presented by Hirsch M.W. 
and Smale S. [4 ]. 

In this chapter, we will establish the condition of the 
occurrence of local bifurcation and Hopf 
bifurcation around each of the equilibrium point of 
a mathematical model proposed by AzharM,and 
Noor Hassan [ 5 ]. 


II. Model formulation [ 5 ]: 

Anecological mathematical model 
consisting of prey-predator model involving prey 
refuge with two different function response is 
proposed and analyzed in [ 5 ]. 
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with initial conditions Wj( 0 ) > 0 , i = 1,2,3,4 . 
Not that the above proposed model has twelve 
parameters in all which make the analysis difficult. 
So in order to simplify the system, the number of 
parameters is reduced by using the following 
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Then the non-dimensional from of system (1) can be 
written as: 


dw 4 
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= w, 


w 2 (l - w 2 ) 


■-a. 


w 4 

= w 1 f 1 (w 1 ,w 2 ,w 3 ,w 4 ) 
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dw 2 

dt 


= w 2 


a 1 w 1 


w 2 


- a 2 - (1 - m)w 3 


= w 2 f 2 (w 4 ,w 2 ,w 3 ,w 4 ) 
dw 3 


( 2 ) 


dt 

dw 4 


= w 3 


a 5 w 4 


= W 4 


a 4 (l — m)w 2 — a 3 — ■ 

a 6 + w 3 . 

= w 3 f 3 (w 4 ,w 2 ,w 3 ,w 4 ) 
a 7 w 3 


.., ' — a 8 = w 4 f 4 (w 1 ,w 2 ,w 3 ,w 4 ) 

dt La 6 + w 3 

with W] ( 0 ) > 0 , w 2 ( 0 ) > 0,w 3 (0) > 

0 and w 4 ( 0 ) > 0 . It is observed that the number 
of parameters have been reduced from twelve 
in the system (1) to nine in the system (2 ) . 
Obviously the interaction functions of the system 
(2 ) are continuous and have continuous partial 
derivatives on the following positive four 

dimensional space. 

| (w 4 ,w 2 , w 3 , w 4 ) 6 R 4 : w 4 ( 0 ) > 0 , 

R+ = |w 2 (0)>0,w 3 (0)>0w 4 (0)>0. 


Therefore these functions are Lipschitzian on R+ , 
and hence the solution of the system (2 )exists and is 
unique. Further, all the solutions of system (2 ) with 
non-negative initial conditions are uniformly 
bounded as shown in the following theorem. 
Theorem 1 : All the solutions of system (2 ) which 
initiate in R' 1 ,.. are uniformly bounded. 


III. The stability Analysis of Equilibrium 
Points of system (2) [5] 

It is observed that, system (2) has at most 
four biological feasible equilibrium point which are 
mentioned with their existence condition in [5 ] as in 
the following: 


l.the Vanishing Equilibrium point: E 0 = 

(0,0,0,0) always exists and E 0 is locally 
asymptotically stable in Int R 4 if the following 
condition hold 

a 2 > 1 (3. a) 

However, it is (a saddle point) unstable otherwise. 
More details see [ 5 ] 


2. The freepredators equilibrium point E x = 

(iVi, m> 2 , 0,0 )exists uniquely in 

Int. /?+ if and only if the following condition 

hold: 


w 2 < ■ 


a 3 


(3.6) 


a 4 (l - m ) 

The following second order polynomial equation 


+ AfA + A 2 — 0 , 

which gives the other two eigenvalues of J (Eg) by: 



Eq is locally asymptotically stable in the R+ . 
However, it is a saddle point otherwise . 


3. The freetoppredator equilibrium pointE 2 = 

(w 1; w 2 , w 3 ,0 ) exists uniquely in the Int. R\if the 
following condition are holds: 

a 3 < min{ a 4 (l — m) , a 4 (1 — m)( 1 

-a 2 )}. (3d) 

The following three order polynomial equation 

[2.3 + R X A^ + R 2 A + R 3 ] ^— cLq + 

= 0 

where 


a 6 + w 3 


--A 


E\ — —( b n + b 22 ) > 0 , 

R 2 — ~(.bnb 22 + b 23 b 32 + b 12 b 2 i), 

R 3 = bnb 23 b 32 > 0 , 

So, either—a R + avW3 — A = 0 . 

a 6 + w 3 


Or2 3 + R 3 A^ + R 2 A + R 3 — 0 

Hence from equation ( 2.9 c ) we obtain that: 

a 7 w 3 

2-2W3 — —a 8 + ~ ’ 

a 6 + w 3 

which is negative if inaddition of condition (2. 5 
h)in [ ] the following condition hold: 


w 3 < 


a 6 a 8 
CZy CLq 


( 3 e ) 


On the other hand by using Routh-Hawirtiz criterion 
equation in [5] has roots ( eigenvalues ) with 
negative real parts if and only if 
R 3 > 0 ,R 3 > 0 and A = R 3 R 2 — R 3 > 0 
Straightforword computation shows that A > 0 
a 3 1 

provided that— -- > - (3/) 

F a 4 (l — ?n) 2 v 

Now it is easy to verify that R x > 0 and/? 3 > 

0 under condition( 2.4 g ).Then all the eigenvalues 

A w ,A W2 and A W3 of equation( 2.9 d )have negative 


real parts. So, E 2 is locally 

4. Finally, the Positive (Coexistence)Equilibrium 

point: 

E 3 = (w 3 ,w 2 ,w 3 , w 4 )exists and it is locally 
asymptotically stable, as shown in [ ]. 


IV. The local bifurcation analysis of 

system( 2 ) 

In this section, the effect of varying the 
parameter values on the dynamical behavior of the 
system (2) around each equilibrium points is 
studied. Recall that the existence of non hyperbolic 
equilibrium point of system ( 2 ) is the necessary but 
not sufficient condition for bifurcation to occur. 
Therefore, in the following theorems an application 
to the Sotomayor’s theorem for local bifurcation is 
appropriate. 

Now, according to Jacobian matrix of 
system( 2 )given in equation( 2.7 ), it is clear to 
verify that for any nonzero vector V = 

(v 4 ,v 2 , v 3 ,v 4 ) t we have: 
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D 2 f(v ,v) = 
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In the following theorems the local 
bifurcation conditions near the equilibrium point are 
established. 

Theorem (3.1): If the parameter a 2 passes through 
the value a 2 — 1 , then the vanishing equilibrium 
point Eq transforms into non-hyperbolic equilibrium 
point and system (2.2) possesses a transcritical 
bifurcation but no saddle-node bifurcation, nor 
pitchfork bifurcation can occur at E 0 . 

Proof: According to the Jacobian matrixJ(E 0 )given 
by Eq. ( 2.8 a ) the system ( 2 ) at the equilibrium 
point E 0 has zero eigenvalue (say A 0x = 0) 

ata 2 — a 2 , and the Jacobian matrixj 0 with a 2 = a 2 
becomes: 

/-a 1 10 0 

Jo ~ J ( a 2 ~ a f> — ( Ul Q 

\0 0 

Now, let v[°l = , p 2 0 ^, vl°\ be the 

eigenvector corresponding to the eigenvalueA 0x — 0 
.Thus (/d ~ A 0 wi l)V^ = 0 , which gives: p 2 ° ] = 

a 4 p| 0 ', v\ 0 ^ = 0 , p 4 °^ = 0 and v[any nonzero real 
number. 

Pet V 7 ^ = ,xp 2 ^ ,^ 3 °' be the 

eigenvector associated with the 

eigenvalue A 0wi = 0 of the matrix JJf. Then 
we have, {jjf — /Iq^/) 1 ! 710 ] — 0 .By solving this 
equation for V 7 ^ we obtain, 

,ip[ ^,0,o) , where xp[ 0 ^ any nonzero 

real number. 

Now, consider: 

df_ = (dfydf^tfk df\ T 

da 2 ° 2 ,CL2 \da 2 ’ da 2 ’ da 2 ' da 2 ) 

= (o, -w 2 , o, o y . 

So, f a 2 (E 0 ,a 2 ) = (0,0,0,0) r and 

hence{jV [ 0 ] ) T f a fE 0 ,a 2 ) = 0 . 

Therefore, according to Sotomayor’s theorem the 
saddle-node bifurcation cannot occur. While the 
first condition of transcritical bifurcation is 
satisfied. Now, since 



Df a? (W ,a 2 ) = 


where Df a2 (W ,a 2 ) represents the derivative of 
fa 2 (W, a 2 )with respect to W = (w 4 ,w 2 ,w 3 , w 4 ) r 
. Further,it is observed that 


Df a2 (E Q ,a 2 )VW = 



O [0] ) [Dfa 2 (E 0 ,a* 2 )vm] = - ai v[ 0] xp[ 0] 

=£ 0 . 

Now, by substitutingin (3.1) we get: 

-2 a^pf'j ^ 

D 2 f(E Q ,a* 2 ){V^ ,V^) = 




/ 


Hence, it is obtain that: 

(y[°]) r D 2 f(E 0 ,a* 2 )(vW ,0°]) 

= -2 a 4 t/;{ 0 ^p| 0 ^ ^ 0 . 

Thus, according to Sotomayor’s theorem (1.15) 
system ( 2 ) has transcrirtical bifurcation atE 0 with 
the parametera 2 = a 2 . 

Theorem ;If the parameter a 3 passes through the 
value a) = a 4 (l — m)w 2 , then the free 
predatorsequilibrium point E 1 = (w, , tv 2 ,0 ,0 ) 
transforms into non-hyperbolic equilibrium point 
and system (2) possesses a transcritical bifurcation 
but no saddle-node bifurcation, nor pitchfork 
bifurcation can occur at E 3 = (w, , w 2 ,0 ,0 ). 

Proof: According to the Jacobian matrix]igiven by 
Eq. ( 2.9 a ) the system ( 2 ) at the equilibrium 
pointEfias zero eigenvalue (sayA 1W3 — 0) at 
a 3 = a 3 , and the Jacobian matrixwith a 3 = a 3 
becomes: 

Ji — J( a 3 — a d) 

/—aq 1 — 2 w 2 0 0 

i ~ a 2 —(1 — m)w 2 — (1 — m)w 2 

0 0 0 0 
0 0 0 —a 6 + a 7 (l - m)w 2/ 

Let l 7 ^ 1 ! = (p-f 1 ^ ,p 2 ^ ,p]^ jP^ 1 ^ be the 
eigenvector corresponding to the eigenvalueA 1W3 = 
0. Thusfa - A 1W3 I)VM =0 
,which gives: 

,,[i] _ C 1 ~ 2a 2 )(l ~ m) [I] tl] 

v 3 > v 2 


CLi 


= — (1 — m)p. 
= 0 , 

where p|^ any nonzero real number. 


3 ^ and p 4 ^ 
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T 

Let ,\p)p be the 

eigenvector associated with the eigenvalueA 1W3 = 
0 of the matrix ]{ T . Then we have, 
(jf — /l llV3 /)¥ / * 1 l = 0 . By solving this equation 

for we obtain, ¥ / l 1 ' = (o , 0 , o) , 

where xp^ any nonzero real number. Now, 
consider: 

T 


3 f_ =f , , = (df L df L d^dU_\ 

da 3 ° 3 ,Cl3 \da 3 ’da 3 'da-i’da*.J 

— (0,0, -w 3 ,0) 7 


and 


So, = (0,0,0,0) r 

hence{f¥^~) f a3 (E 1 ,a 3 ) = 0 . 

Therefore, according to Sotomayor’s theorem (1.15) 
the saddle-node bifurcation cannot occur. While the 
first condition of transcritical bifurcation is 
satisfied. Now, since 

/O 0 0 0 

n f (va ~ a / 0 0 0 0 

D fa 3 (W,a 3 ) - I 0 0-10 

Vo 0 0 0' 

where Df a3 (W ,a 3 ) represents the derivative of 
f a ( W , a 3 )with respect to W = (w 1 ,w 2 ,w 3 , w 4 ) T 
.Further,it is observed that 


Df a (E 4 ,af)V^ = 


/° 

0 

0 

0 \ 

7i [11 \ 

v\ 1] 

[ 0 

0 

0 

0 1 

1 ° 

\o 

0 

0 

-1 

0 

7 

0 / 

r> 1] 

w 


0 

0 

—v. 


[i] 


0 


V 


0 [1] ) [Df a3 (E 1 ,a* 3 )V^]=-vf ] ff 
* 0 . 

Now, by substituting]/^ in (3.1) we get: 
D 2 f(E 4 ,af){V^ ,V^) 

( 2 (1 - m.y(v\ 1] ^ ^ 

2 (1 — m) 2 ^ 1 ^ 

-2 a 4 ( 1 - m) 2 ^ 1 ^) 

0 

Hence, it is obtain that: 

(ipWyD 2 f(E 1 ,a* 3 )(V ll l ,V ll] ) = 

-2 a 4 (l - m) 2 ^ 0 . 

Thus, according to Sotomayor’s theorem (1.15) 
system (2) has transcrirtical bifurcation atE^ with 
the parametera 3 = a 3 

Theorem (3.3 ):If the parameter a 8 passes through 




the. value a 3 = 


a-j W 4 

a 6 +w 3 


, then the free top predators 


equilibrium point E 2 = (wi , w 2 , w 3 ,0 ) transforms 
into non-hyperbolic equilibrium point and system 
(2.2) possesses a transcritical bifurcation but no 


saddle-node bifurcation, nor pitchfork bifurcation 
can occur at E 2 = (w, ,w 2 , w 3 ,0 ). 

Proof: According to the Jacobian matrixj 2 given by 
Eq. ( 2.10 a ) the system ( 2.2 ) at the equilibrium 
pointE 2 has zero eigenvalue (say A 2w4 = 0) at 

a 8 — a g , and the Jacobian matrix J 2 with 
a 8 — a 8 becomes: 

1*2 =/(«8 = aS) = [ e ij] 4x4 where 

n n = —a 1 ,n 12 =1-2 w 2 ,n 13 = 0,n 14 = 0,n 21 

—— flj , 

n 22 = ~a 2 - (1 - m)w 3 ,n 23 = -(1 - m)w 2 ,n 24 
= 0, 

n 3 i = 0 ,n 32 = a 4 (l - m)w 3 ,n 33 

= —a 3 + a 4 ( 1 - m)w 2 , 

a 5 w 3 

n 34 = -rrr" ,n 41 = 0,n 42 = 0,n 43 = 0,n 44 

a 6 + w 3 

= 0 . 

Let = (v[ 2 \v^ ,v^ be the 

eigenvector corresponding to the eigenvalueX 2v/3 = 
0 .Thus(j 2 — 3 2w3 l)V^ 2 ^ — 0 ,which gives: 

[ 2 ] __ a 5 [a 4 (l — m) — 2a 3 ] _ j 2 ] 


a x a 6 a 4 (l — m) 2 (l — a 2 ) — a 4 a 3 a 4 


J21 


a 5 (l-m) 


and v. 


[ 2 ] 


a 4 a 6 (l - m) + [a 4 (l - m) - a 3 ] 


,[ 2 ] 


a 5 


J2] 


a 4 a 6 (l - m) + [a 4 (l - m) — a 3 ] 4 
wherev 4 2 ^ any nonzero real number. 

Let = (fpy ,tpy ,tpy be the 

eigenvector associated with the 

eigenvalue 3 2w3 = 0 of the matrix J 2 T .Then we 
have, (j 2 T — 3 2w J ) ’f'* 2 ] = 0 . By solving this 

equation forf^we obtain, ¥^ 2 ' = (o,0,0,t/7) , 
[21 

where xp 4 1 any nonzero real number. Now, 
consider: 

df_ =f , , = (dA fh dh_ df\ T 

da 8 a ® ’ as \da 8 ' da 8 ’ da 8 ’ da 8 ) 

= (0,0,0, -w 4 ) T . 

So, fa 8 (V 2 > a fi) — (0,0,0, 0) T and 

hence{V [2] ) T f ae (E 2 ,a 8 ) = 0 . 

Therefore, according to Sotomayor’s theorem (1.15) 
the sacldle-node bifurcation cannot occur. While the 
first condition of transcritical bifurcation is 
satisfied. Now, since 

20 0 0 ON 

0 0 0 0 
0 0 0 0 
^0 0 0 -V 

where Df ag (w ,a 8 ) represents the derivative of 

f ag ( X, a 8 )with respect to w = (w 4 , w 2 , w 3 , w 4 ) T 
. Further, it is observed that 


Df a yx,a 8 ) = 
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Df a JE 2 ,a* 8 )VW = 


/° 

0 

0 

0 \ 

/^ [21 \ 

7 i 21 

1 0 

0 

0 

0 1 

1 ° 

Vo 

0 

0 

0 

0 

-V 

v? ] 

U/ 


0 

0 

0 

k-v, 


[ 2 ] 


0 [2] )' [Df as {E 2 ,al)V p]] = -1 4 2 V [2] 


* 0 


Now, by substitutingin (3.1) we get: 
D 2 f(E 2 ,CL* 8 )(yW,VW) 

I A’?') 1 


2v. 


[ 2 ] 


—2(1 — m)v^vlf 
a 4 ( 1 — m)v P' 


[ 2 ] 


+ 


&6 ^7 ^4 


(a 6 + w 3 ) 2 


V 


</4 2] * 0 


[ 2 ] [ 2 ] 
a b a 7 v\ v\ 

(a 6 + w 3 ) 2 

Hence, it is obtain that: 

(yW) T D 2 f(E 2 ,a* 8 )(yw 

_ o a 6 a 7 ^ 2] ^ 2] 

(06 + W 3 ) 2 ^ 

Thus, according to Sotomayor’s theorem (1.15) 
system (2.2) has transcrirtical bifurcation 
atE 2 with the parametera 8 = a 8 . 

Theorem (3.5): Suppose that the following 
conditions 

a 7 > a 8 (1 + (1 - m)a 6 ) (3.4) 

are satisfied. Then for the parameter value 

a 2 = 1 — (1 — m) - 

CLj CIq 

system (2.2) at the equilibrium point E 3 — 
(w( ,w 2 ,w 3 ,w|) has saddle-node bifurcation, but 
transcrirtical bifurcation, nor pitchfork bifurcation 
can occur at E 3 . 

Proof: The characteristic equation given by 
Eq.( 2.12 b )having zero eigenvalue ( say A 4 — 0 ) 
if and only if C 4 = 0 and then E 3 becomes a 
nonhyperbolic equilibrium point. Clearly the 
Jacobian matrix of system( 2.2 )at the equilibrium 
point E 4 with parameter a 2 = a 2 becomes: ]\ = 
J(a 2 = af) = [£y] 4x4 where 
£ ij = d-ij for all i,j = 1,2,3,4 excepts^ which is 
given by: 

CLc. CLj 

£ 22 = -a* - (1 ~m) - . 

(Zy CLq 

Note that, af > 0 provided that condition (3.4) 
holds. 

Let j/Pl = ,v^ be the 

eigenvector corresponding to the eigenvalueA 4 = 
0 . ThusQl — A 4 I) = 0, which gives: 


P] 

V, — 


,[3] _ 


1-2 w 2 * 


-V- 


[3] „[3] 


V, = 


(d 12 + d 22 ) 


22' [3] 

, 


a 4 ‘ J (1 — m)wf “ 
~(^i 2 + d 22 )a e a 7 wf _ j 3 ] 


d 23 (a 6 + w 3 *)(a 6 a 8 + (a 8 - a 7 )w 3 *) 
andv 7 ' any nonzero real number. Herevf, aj 3 
not equal zero under the existence condition 

(2.5t). 

Let ,tp 2 \tp^ be the 

eigenvector associated with the 

eigenvalueA 4 — 0 of the matrix])) .Then we have, 
(/f — 2 4 /) l / y3 = 0. By solving this equation 
foiW^ 3] we obtain: 

I//[3] 


/ 


( d 


12 


v 4 3] 

</4 3] 

+ d 22 ) 




[3] 


a 4 (l - m)w 3 

1_ ( rf 12 + d 22 ) _ [3] 

\a 4 (l - m)(a 6 a 8 + (a 8 - a 7 )w 3 *) 2 / 
[4] 


So, 


where xp 2 1 any nonzero real number. 

Now, 

*L = t 32 (M/,a 2 ) = f^,^,^,^) r = 

da 2 Ja2 ^ \3a2 9a2 3a2 3a2/ 

(0 , -w 2 ,0,0) r . 

)a 2 (£ 3 , af) = (0 , —w 2 ,0,0) r and hence 

(^ [ 3])7 a2 (O,a 2 0 = -w|^ 41 
=£ 0 . 

Therefore, according to Sotomayor’s theorem (1.15) 
the saddle-node bifurcation occure at E 3 under the 
existence condition (2.5 q ). 

The Hopf bifurcation analysis of system ( 2 ) 

In this section, the possibility of occurrence 
of a Hopf bifurcation near the postiveequilibrium 
point of the system ( 2 ) is investigated as shown in 
the following theorems. 

Theorem (3.8): Suppose that the following 
conditions are satisfied: 

max f—Z 6 , (i3+rf23ti32) | < a 4 < min F where 

t ^33 J 

F= 

Id 


h -£3 h Q5 a 6 a 7 u '3 H '4 ^33(~ a l ti 12+^33) ) 

33 ’3(l-m)W3 ’ d 34 d 43 ’^(ag+wf) 3 ’ d33+(3 J 

max [l 3 , a 4 l 3 ,2d 33 Z 4 + 3d n } < < 


Co < min 


CiC 2 C 2 


u 


(3.9) 


-3 

C 1 3 -4A 1 >0 (3.10) 

Then at the parameter value a 2 — a 2 , the 
system ( 2.2 ) has a Hopf bifurcation near the 
pointE 3 . 

Proof: According to the Hopf bifurcation theorem 
(1.17), for n=4, the Hopf bifurcation can occur 
provided thatC i (a 2 ) > 0 ; 
i = 1,3 , A 4 > 0 , 
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C( — 4 A,> 0 andA 2 (d 2 ) = 0. Sti'aight fonvard 
computation gives that: if the following 

conditions (2.10 d), (2.10 e) ,(3.7),(3.8),(3.9) and 
(3.10) are hold, thenCfaf) >0; i = 1,3, A,> 
0 and Cj 3 — 4 A 1 > 0 . 

Now, to verify the necessary and sufficient 
conditions for a Hopf bifurcation to occur we need 
to find a parameter satisfy A 2 = 0. Therefore, it is 
observed thatA^= 0 gives that: 

C 3 (C 1 C 2 - C 3 ) - Cl C 4 = 0 
Now, by using Descartes rule Eq.(3.11) has a 
unique positive root saya 2 = d 2 such that: 

£ 1 a 2 3 +X 2 “ 2 2 +L 3 a 2 +£ 4 = 0 


A 4 y — OCy + t(X 2 ) ^-4 2 — 0*1 



Clearly, Re(A 4k (a 2 )) I a2=S2 = afa 2 ) = 0 ,k = 
1,2 that means the first condition of the 
necessary and sufficient conditions for Hopf 
bifurcation is satisfied ata 2 — a 2 .Now, according 
to verify the transversality condition we must 
prove that: 

0(a_ 2 )na 2 ) + f(a 2 )0(a 2 )^ 0 , 

whereO , V, TancfP are given in (1.29). Note that 


I'l = hh < 0 

L 2 — (1 — m)w 3 (f — d 33 )Z 2 

— d 34 d 43 (d 34 d 43 + 3d n 
+ 2 ^d 33 ) 

+r 3 Z 2 + dii Cdn f 3 — d 12 d 21 d 33 — d 33 lf) 


for a 2 = a 2 we have a 4 = 0 and a 2 = I— , 

substituting into (1.29) gives the following 
simplifications: 


— d 12 d 21 d 33 Z 1 . 

L 3 = -(1 - m) 2 wfl 4 l 2 

+ (1 

— m)w 3 (d 34 d 43 d n (3 + l t j) 

+2 d 34 d 43 (d 34 d 43 — d 33 l 2 ) — 2 T 3 l 2 — 2 d 41 l 3 

+ d 12 d 21 d 33 l 4 

-r 3 d n (l 3 + l 2 ) + d 12 d 21 d 33 (Z 5 + 3(1 - m)w 3 ) 
+d 33 L ( 2 d 41 — r 3 + 2(1 - m)w 3 d n ) 

£ 4 = (l-m)w 3 T 3 I 1 a 2 +d 12 d 21 ) 

- (1 - mYwffl 2 

+r 3 ^12^21(^11 + ^33) + d 33 d 34 d 43 (d 12 d 21 — d 41 ) 

+ (1 - m)w 3 d 33 (l 2 + Z 4 
— Zid 14 ) 

— (1 - m) 2 w 3 * 2 (d 11 (Z 5 

— d 33 d 12 d 21 — d 34 d 43 ) 

— d 34 d 43 (d 34 d 43 + d 41 

^3^11(^34^43(^1 ~ ^11^3) 

+ d n ) + d 41 r 3 d 23 d 32 

— d 11 d 33 (r 3 i 6 — d 33 d 12 d 21 Z 4 

— ^ 3 ^ 11 ) 

+ (1 - m)w 3 d n d 34 d 43 (Z 3 
+ d 34 d 43 ). 


^(d,) - —2 C 3 (a 2 ) , 

_ a 2 (a 2 ) 

0(d 2 ) = 2 -^-(^^-2^) , 

L i 

Co , 

0 (^ 2 ) — Q( a 2 ) — 7^02 ( a 2 ) < 

L 1 


r(a 2 ) = a 2 (a 2 ) ( C 3 (a 2 ) - — Ci(a 2 ) 


where 



dCy I 


Cl = 

da 2 ° 2= 

a 2 1 ' 

dC 2 

1 


d(X 2 

'a 2 =a 2 ~ 

~fl . 

dC 3 

1 


dci 2 

'a 2 =a 2 

~fl > 


dC 4 


c\ = 

da 2 ° 2 

=CL 2 = d ll 


Then by using Eq. (1.30) we get that: 
0(d 2 )na 2 ) + f(d 2 )0(d 2 ) = 


2C 3 (a 2 ) (^11^34^43 + 
2C3(12+C3C1) . 


c 3 l 2 \ _ 2 a 2( a 2o 
Ci / Ci 


(CiC 2 - 


andly 

= d 14 

+ ^33 

> 0 , 






l 2 = d 

34^43 “ 

- d n 

d 33 > 0 


h = 

^33 ~ 

d 23 d 32 

> 0 , 

U 

= d 34 d 43 

+ d 12 d 21 



< 0 




h = 

^33 — 

d 12 d 2 i 

> 0 , 

h 

= d 12 d 21 

+ d 23 d 32 


< 0 

/Von/, afa 2 = a 2 the characteristic equation can be 
written as: 

( A|+ |)(^ + ClA4+ ^)=° ' 
which has four roots 

*41,2 = ±i^andX 43A = i(-Ci ± Jcf^ . 

Clearly, ata 2 = a 2 there are two pure imaginary 
eigenvalues (A 41 and A 42 ) and two eigenvalues 
which are real and negative. Now for all values 
ofa 2 in the neighborhood ofa 2 , the roots in general 
of the following form: 


Now, according to conditions (3.8) and (3.9) we 
have: 

0(d 2 )nd 2 ) + f(d 2 )0(d 2 ) * 0 . 

So, n/f? obtain that the Hopf bifurcation occurs 
around the equilibrium point E 3 at the 
parametem 2 = a 2 . 

■ 

V. Numerical analysis of system( 2 ) 

In this section, the dynamical behavior of 
system ( 2 ) is studied numerically for different sets 
of parameters and different sets of initial points. The 
objectives of this study are: first investigate the 
effect of varying the value of each parameter on the 
dynamical behavior of system (2) and second 
confirm our obtained analytical results. It is 
observed that, for the following set of hypothetical 
parameters that satisfies stability conditions of the 
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positive equilibrium point, system (2) has a 
globally asymptotically stable positive equilibrium 
point as shown in Fig. ( 1). 

a 3 — 0.7 ,a 2 = 0.2 ,a 3 = 0.1 
, a 4 = 0.5 , a 5 = 0.6, m = 0.8 ( 14 ) 
a 6 = 0.4 ,a 7 = 0.5 ,a 8 = 0.2 
,m = 0.8 (14) 




1 




Fig. ( 1 ): Time series of the solution of 
system ( 2 )that started from four different 

initialpoints(0.3,0.4,0.5,0.6), (2.5,0.5,2.5,1.5), 
(2,1.5,2.5,2.5) and (2.5,2.5,1.5,0.5) for the data 
given by( 2.1). {a) trajectories ofw 1 as a function 
of time, (b) trajectories ofw 2 as a function of 
time, (c) trajectories ofw 3 as a function of time, 
(d) trajectories of w 4 as a function of time . 

Clearly, Fig. (1 ) shows that system (2) has a 
globally asymptotically stable as the solution of 
system ( 2 )approaches asymptotically to the positive 
equilibrium point E ? = ( 0.34,0.6,0.4,0.08 ) 


starting from four different initial points and this is 
confirming our obtained analytical results. 

Now, in order to discuss the effect of the parameters 
values of systemf 2) on the dynamical behavior of 
the system, the system is solved numerically for the 
data given in( 14 ) with varying one parameter at 
each time. It is obser\>ed that for the data 
given in( 14 )with 0.1 < a t < 1, the solution of 
system( 2 )approaches asymptotically to the positive 
equilibrium point as shown in Fig.(2)for typical 
value a 1 = 0.2 . 



Fig. (2 ): Time series of the solution of system (2) 
for the data given by{14) with a } = 0.2, which 
approaches to {1.2 ,0.6 ,0.4,0.08) in the interior 
of <. 

By varying the parameter a 2 and keeping the rest 
of parameters values as in (14) , it is observed 
that for 0.1 < a 2 < 0.44 the solution of system 
(2) approaches asymptotically to a positive 
equilibrium point E 3 , while for 0.44 < a 2 < 
0.6 the solution of system (2) approaches 
asymptotically to E 2 = (w, ,w 2 ,w 3 ,0) in the 
interior of the positive quadrant ofw 1 w 2 w 3 — 
plane as shownin Fig. (3) for typical 

value a 2 = 0.5 . 



Fig.( 3 ): Times series of the solution of system 
( 2 )for the data given by ( 14 )with a 2 = 0.5 . 
which approaches to( 0.34,0.4,0.2,0 )in the 
interior of the positive quadrant of w, w 2 w 3 — 
plane. 

while fori). 6 < a 2 < 1 the solution of system 
(2) approaches asymptotically to E 3 = 
(w, , w 2 ,0,0 ) in the interior of the positive 
quadrant ofw 1 w 2 —plane as shown in Fig. 
(4 )for typical value a 2 = 0.8. 



Fig.(4):Times series of the solution of system (2) 
for the data given by (14) with a 2 — 0.8 which 
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approaches to( 0.23,0.2,0,0)in the interior of the 
positive quadrant of w l w 2 — plane . 

On the other hand varying the parameter a 3 and 
keeping the rest of parameters values as in( 14 ), 
it is observed that for 0.1 < a 3 < 0.145 the 
solution of system^ 2 ) approaches asymptotically 
to the positive equilibrium point E 3 . while 
/b/0.145 < a 3 < 0.2 the solution of system ( 2 ) 
approaches asymptotically to 

E 2 = (u>i , w 2 , w 3 ,0 ) in the interior of the 
positive quadrant of w, w 2 w 3 — plane, for 
0.2 < a 3 < 1 the solution of system ( 2 ) 
approaches asymptotically to E 1 = (w, , w 2 -0,0 ) 
in the interior of the positive quadrant of 
w x w 2 — plane. 

Moreover, varying the parameter a 4 and keeping 
the rest of parameters values as in ( 14 ), it is 
observed that for 0.1 < a 4 < 0.26 the solution of 
system (2) approaches asymptotically 
toE 1 = (w, , w 2 ,0,0 ) in the interior of the 
positive quadrant of w 1 w 2 — plane, while 
for 0.26 < a 4 < 0.4 the solution of system ( 2 ) 
approaches asymptotically to 

E 2 — (w, , w 2 , w 3 ,0 ) in the interior of the 
positive quadrant ofw A w 2 w 3 — plan.and for the 
0.4 < a 4 < 1 the solution of 

system ( 2 ) approaches asymptotically to the 
positive equilibrium point E 3 . 

For the parameter 0.5 < a 5 < 1 the solution of 
system (2) approaches asymptotically to a 
positive equilibrium point E 3 . 

For the parameter 0.1 < a 6 < 0.2 the solution of 
system (2) approaches asymptotically to 
E 2 = (w, , w 2 , w 3 ,0 ) in the interior of the 
positive quadrant ofw 1 w 2 w 3 — plan. While for 
the 0.2 < a 6 < 0.5 the solution of system ( 2 ) 
approaches asymptotically to a positive 
equilibrium point E 3 . 

For the parameters values given in( 14 ) with 
0.4 < a 7 < 0.6 the solution of system ( 2 ) 
approaches asymptotically to the positive 
equilibrium pointE 3 . 

For the parameter 0.1 < a 8 < 0.32 the solution 
of system ( 2 ) approaches asymptotically to the 
positive equilibrium point E 3 . while for the 
0.32 < a 8 < 1 the solution of system ( 2 ) 
approaches asymptotically to 

E 2 = (w, , w 2 , w 3 ,0 ) in the interior of the 
positive quadrant of w, w 2 w 3 — plan. 

Moreover, varying the parameterm and keeping 
the rest of parameters values as in( 14 ), it is 
observed that for 0.1 <m< 0.7 the solution of 
system ( 2 ) approaches asymptotically to the 
positive equilibrium point E 3 . while for 
0.7 < m < 0.75 the solution of system (2 ) 
approaches asymptotically to 

E 2 = (w, , w 2 , w 3 ,0 ) in the interior of the 
positive quadrant of w 1 w 2 w 3 — plane, and 
foiO.75 < m < 1 


the solution of system ( 2 ) approaches 
asymptotically to E x = (w, , w 2 ,0,0 ) in the 
interior of the positive quadrant of w, w 2 — 
plane. 

Finally, the dynamical behavior at the vanishing 
equilibrium point E () = (0,0,0,0) is 
investigated by choosing a 2 = 2 and keeping 
other parameters fixed as given in ( 14 ), and 
then the solution of system ( 2 ) is drawn in Fig. 
(5). 



FigJ 5 ): Times series of the solution of system ( 2 ) 
for the data given by (14) of tut h a 2 = 2, 
which approaches ( 0, 0, 0, 0). 

VI. Conclusions and discussion: 

In this chapter, we proposed and analyzed 
an ecological model that described the dynamical 
behavior of the food chain real system. The model 
included four non-linear autonomous differential 
equations that describe the dynamics of four 
different population, namely first immature prey 
mature prey (iV 2 ), mid-predator (iV 3 ) and 
(1V 4 ) which is represent the top predator. The 
boundedness of system (2.2) has been discussed. 
The existence condition of all possible equilibrium 
points are obtain. The local as well as global 
stability analyses of these points are carried out. 
Finally, numerical simulation is used to specific the 
control set of parameters that affect the dynamics of 
the system and confirm our obtained analytical 
results. Therefore system (2.2) has been solved 
numerically for different sets of initial points and 
different sets of parameters starting with the 
hypothertical set of data given by Eq.( 4.1 ) and the 
following observations are obtained. 

1. System (2 ) has only one type of attractor in 
Int. R\ approaches to globally stable point. 

2. For the set hypothetical parameters value given in 

(14 ), the system (2) approaches asymptotically to 
globally stable positive point E 3 = 

( 0.34,0.6,0.4,0.08 ).Further, with varying one 
parameter each time, it is observed that varying 
the parameter values, a t i = 1,5 and 7 do not 
have any effect on the dynamical behavior of 
system (2) and the solution of the system still 
approaches to positive equilibrium pointy = 
(w 4 , W 2 , W 3 , w|). 

3. As the natural death rate of mature prey a 2 
increasing to 0.43 keeping the rest of parameters 
as in eq.(14),the solution of system (2) 
approaches to positive equilibrium point f; 3 . 
However if 0.44 < a 2 < 0.6 , then the top 
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predator will face extinction then the trajectory 
transferred from positive equilibrium point to 
the equilibrium point E z = 

(wj ,w 2 ,w 3 ,0) moreover, increasing a 2 > 0.6 
will causes extinction in the mid-predator and 
top predator then the trajectory transferred from 
equilibrium point E 2 = (w 4 , w 2 ,w 3 ,0 ) to the 
equilibrium point E 4 = (w 4 , w 2 ,0,0 ), thus, thea 2 
parameter is bifurcation point. 

3. As the natural death rate of top predator a 3 

increasing to 0.144 keeping the rest of 
parameters as in eq.(14),the solution of system 
(2) approaches to positive equilibrium pointy. 
However if 0.145 < a 3 < 0.2, then the 

top predator will face extinction then the 
trajectory transferred from positive equilibrium 
point to the equilibrium point E 2 = 
(w 4 , w 2 ,w 3 ,0 ) moreover, increasing 0.2 < a 3 
will causes extinction in the mid-predator and 
top predator then the trajectory transferred from 
equilibrium point E 2 = (w 4 , w 2 ,w 3 ,0 ) to the 
equilibrium point E x = (m\ , w 2 ,0,0 ), thus, thea 3 
parameter is bifurcation point. 

4. As the natural death rate of top predator a 4 

increasing to 0.1 keeping the rest of parameters 
as in eq.(14),the solution of system (2) 
approaches to E x — (w 4 , w 2 ,0,0 ). However if 
0.26 < a 4 < 0.4, then the trajectory transferred 
from E 3 = (w 4 , w 2 ,0,0 ) to the equilibrium 
point E 2 = (w 4 , w 2 ,w 3 ,0 ) moreover, 

increasing a 4 > 0.4 will then the trajectory 
transferred from equilibrium point E 2 = 
(w 4 , w 2 ,w 3 ,0 ) to the positive equilibrium point 
E 2 — (w 4 , w 2 , , w 4 ), thus, the a 4 parameter is 

bifurcation point. 

5. As the half saturation rate of top predator a 6 

decreases under specific value keeping the rest 
parameters as in eq( 14 ) the top predator will 
face extinction and the solution of system (2) 
approaches to E 2 = (w 4 ,w 2 ,w 3 ,0). but where 
a 6 increases above specific value the trajectory 
transferred from E 2 = (w 4 , w 2 , w 3 ,0 ) to the 
positive equilibrium point E 3 = 
(w 4 ,w 2 , , w 4 ),Thus, the parameter a 6 is a 

bifurcation point. 

6. As the natural death rateof the top 
predatora 8 >0.1 increasing keeping the rest of 
parameters as in eq( 14 ) system (2) has 
asymptotically stable positive point in 
Int.f? 4 .However increasing a 8 > 0.32 will causes 
extinction in the top predator then the trajectory 
transferred from positive equilibrium point to 
the E 2 = (w 4 ,w 2 , w 3 ,0 ). Thus, the parameter a s 
is a bifurcation point. 

7. As the number of prey inside the refuge m 

increasing to 0.68 keeping the rest of parameters 
as in eq.(14), the solution of system (2) 
approaches to positive equilibrium point E 3 . 
However if 0.7 < m < 0.75 then the 


top predator will face extinction then the 
trajectory transferred from positive equilibrium 
point to the equilibrium point E 2 = 

(w 4 , w 2 ,w 3 ,0 ) moreover, increasing 0.75 < m 
will causes extinction in the mid-predator and 
top predator then the trajectory transferred from 
equilibrium point E 2 = (w 4 , w 2 ,w 3 ,0 ) to the 
equilibrium point E 4 = (w 4 , vv 2 ,0,0 ), thus, them 
parameter is bifurcation point. 
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